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Abstract 

In this work, analytical solutions were derived (for any methanol oxidation reaction order) for the 
profiles of methanol concentration and proton current density by assuming diffusion mass transport 
mechanism, Tafel kinetics, and fast proton transport in the anodic catalyst layer of a direct methanol 
fuel cell. An expression for the Thiele modulus that allows to express the anodic overpotential as a 
function of the cell current, and kinetic and mass transfer parameters was obtained. For high cell current 
densities, it was found that the Thiele modulus {(fy^) varies quadratically with cell current density; 
yielding a simple correlation between anodic overpotential and cell current density. Analytical solutions 
were derived for the profiles of both local methanol concentration in the catalyst layer and local anodic 
current density in the catalyst layer. Under the assumptions of the model presented here, in general, 
the local methanol concentration in the catalyst layer cannot be expressed as an explicit function of the 
position in the layer. In spite of this, the equations presented here for the anodic overpotential allow the 
derivation of new semi-empirical equations. Keywords: Electrochemical kinetics; Catalyst layer; Mass 
transport; Methanol electro-oxidation; Proton transport. 

1 Introduction 

The direct methanol fuel eell (DMFC) has been proposed as an alternative clean and safe power source, due to 
the easy transportation and storing of liquid methanol and the corresponding low carbon dioxide emissions of 
the cell. However, the DMFC has performance inconveniences attributed to the methanol crossover through 
the solid polymer electrolyte (SPE) and kinetic limitations as the poisoning caused by species like carbon 
monoxide at the anode. 

The simplest approach to the DMFC polarization curve is by using semi-empirical equations, and along with 
the nonlinear least-squares minimization technique, the polarization curve of a fuel cell (including a DMFC) 
can be fitted. One of these semi-empirical equations was first proposed by Srinivasan et. al [1]: 

E = Eo-b\n I - RI (1) 

Eo=E, + bliiIo (2) 

where / is the cell current density and E is the cell potential, and the other quantities are adjustable 
parameters. The equation of Srinivasan et al. [1] is not suitable when mass transfer limitations occur. On 
the other hand, Kim et al. [2] showed that the equation 

E ^ Eo-blnl - RI -m exp(n/) (3) 

accurately fits the polarization curve in the mass transfer region where the equation of Srinivasan et al. [1] 
is not applicable; however, Eq. (3) fails in the low cell current region. Considering this, Squadrito et al. [3] 
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proposed a semi-empirical equation in which the term mexp(n/) in Eq. (3) was replaced by aJ*^ ln(l — (51); 
this modified equation works well over the whole cell current range. Additionally, semi-empirical models 
for low temperature fuel cells that are based on mass transport in the DMFC have been published. For 
example, by invoking an empirical mass transfer coefficient and Tafcl kinetics, Argyropoulos et al. [4] 
derived an equation with adjustable parameters. This equation accurately fits the polarization curve even 
in the low current region where other theoretical models fail; however, the semi-empirical equation derived 
by Argyropoulos et al. [4] does not take into account the effects of polarization caused by the reduction 
of oxygen in the cathode. For this reason, Tu et al. [5] proposed a semi-empirical equation that allows to 
include the effects in the polarization curve of both anodic and cathodic ovcrpotcntials, and the effect caused 
by the mixed potential in the cathode. 

Equations derived from mathematical modeling are more rigorous than semi-empirical equations and yield 
more reliable results. One of the objectives of mathematical modeling is the prediction of fuel cell polarization 
curve, in which the effects of activation, and Ohmic and concentration overpotentials can be observed. 
Estimation of the effects of parameters is another objective of mathematical modeling. For this reason, it is 
necessary to propose models and perform simulations that make it possible to establish such effects on the 
fuel cell performance. 

Baxter et al. [6] considered the anode as a porous electrode covered by a selective SPE. By taking into 
account the transport of methanol, carbon dioxide, water and protons, and Butler- Volmer kinetics, a pseudo 
two-dimensional model was proposed, in which the catalyst layer consists of a liquid layer, a bond layer and 
a matrix layer. Baxter et al. [6] considered the species exchange between the bond layer and the liquid layer; 
this mass transfer process is characterized by effective mass transfer coefficients. The results of Baxter et al. 
[6] suggest that the anode thickness, the surface area, and the charge transfer coefficient are important for 
the determination of the anode behavior. 

Scott et al. [7] developed a model for the liquid-fed DMFC which considers the influence of the carbon 
dioxide, methanol and water transport at the fuel cell anode. The water transport is assumed to be caused 
by electro-osmotic drag and molecular diffusion. Those authors simulated the fuel cell polarization curve 
using a semi-empirical equation for the open circuit potential and a capillary pressure model. In this case, 
the spatial variation of methanol concentration in the catalyst layer was not considered. 
Kulikovsky [8] derived an analytical solution for proton transport across the cathode of a polymer electrolyte 
fuel cell (PEMFC), Tafel kinetics and no variation of the oxygen concentration throughout the cathodic 
catalyst layer were assumed; these solutions are appropriate for low cell current densities. 
Wang et al. [9] proposed a model for oxygen reduction at the electrocatalytic layer, in which the oxygen 
and proton transport occurs through spherical agglomerates, to study the case of fast proton transport, fast 
oxygen diffusion, and mixed control. For the case of low protonic conductivities, the results of Wang et al. 
[9] suggest that the active part of the catalyst layer is located at the region close the membrane. 
Many works have been published on modeling and simulation of fuel cells and are based on numerical methods 
[6, 9, 10, 11, 12] to solve the transport equations, only few of those works have implemented analytical 
methods to simulate the fuel cell. However, the simulation of fuel cells does not dependent on the procedure 
of calculation, and analytical methods do not offer computational advantages over numerical methods. In 
spite of this, procedures based on analytical solutions contribute to the analysis and understanding of the 
processes that occur in the electrodes of fuel cells. 

In this work, we have developed a one-dimensional analytical model of the anode side of a DMFC. Fick's 
law for methanol transport in the anode, and Tafel kinetics for methanol electro-oxidation, and fast proton 
transport were assumed. A non-linear differential diffusion-reaction equation for methanol transport, which 
was solved analytically, was obtained under these assumptions. The solution so obtained is an implicit 
function of the methanol concentration at the catalyst layer and it must be solved by an iterative procedure. 

2 Theoretical model 

In this model it is assumed that: 

1. The fuel cell temperature is constant. 
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2. The pressure difference between tlie botfi sides of tlie cell is negligible. 

3. The membrane is fully hydratcd. 

4. The effects caused by carbon dioxide evaporation, methanol evaporation and water are negligible. 

5. Carbon monoxide poisoning is negligible. 

6. Steady state operation conditions. 

7. Constant cathodic ovcrpotcntial. 

8. Methanol transport throughout the electrolyte phase in the catalyst layer is neglected. 



2.1 Methanol transport in CL 

Figure 1 shows a schematic representation of the anode side, this is divided into four zones; flow channel, 
diffusion layer (DL), catalyst layer (CL), and the membrane. The domain of the z variable is the region 
< z < 5ci in the CL. CL consists of a carbon supported solid catalyst. A finite external load induces three 
basic phenomena in the CL: (i) methanol transport, (ii) proton transport, and (iii) a finite overpotential 
on the anode. Methanol molecules in the DL are attracted towards the CL by electro-osmotic drag, and a 
concentration gradient caused by methanol electro-oxidation, once the molecules reach the CL, the electro- 
oxidation of methanol molecules decrease the methanol flux throughout the CL. In the membrane, methanol 
flux is considered constant. From the molecular diffusion flux deflnition (in steady state) it follows that 

-Vt,^^''^+v{z)Ct[z)^Ntk{z) (4) 

here, C£j is the local methanol concentration at the CL, z; is a convective velocity [14], Tf^"^ is the methanol 
effective diffusion coefficient in the CL and it is given by the Bruggemann relation 

vtf^4^V^ (5) 

The methanol flux is expressed as follows 

iV^(z) = ^^^^+7V^f (6) 

here, / is the cell current density, j is the local anodic current density, is the local methanol flux at 
the CL, N'^ is the methanol crossover flux, and tt-m is the number of electrons transferred per methanol 
molecule. 



2.2 Analytical solutions for methanol transport in the CL 

Given that the convective velocity in the CL is small (no pressure gradients), the transport of methanol by 
convection mechanism in the CL can be neglected, and hence, Eq. (4) can be expressed as follows 

where the methanol flux at the DL: 

< = -^:+^M (8) 
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The volumetric anodic current {dj /dz) is given by 



where 5* is the catalyst specific superficial area, f^^ is the reference current, is the reference methanol 
concentration, a is the charge transfer coefficient, 7 is the methanol electro-oxidation reaction order, and 77 
is the anodic overpotential (assumed independent of z). Methanol flux is related with volumetric current as 
follows: 



d< 1 dj 



dz timJ' dz 

substitution of this equation into Eq. (9) gives 

7 



dz^ nM-F<^^ ^Vi?r Cm 
Define the dimensionless variables 

, ^ r(n-^MiO , j(C) 
I- ^ Q(C) - 

nuTkciC^^'^^ dC 



The Thiele modulus is defined as 



where C^^*^' is the methanol concentration at the CL/membrane interface, and 

-T-*cl,cfF 



Define 

^o = Q(0) = — ^ Oi = '3(i) = — ^ 



and 



here CJJJ is the methanol concentration at the CL/membrane interface. Note that 



(On 



C(0) 



cl/dl 
M 



^cl/dl 



Using the previous dimensionless variables we can write Eq. (9) as 

dJ_(t^ 
'dC ~ J 



and Eq. (11) as 



(10) 



^^^exp(^,) £0) ^0 (11) 



(12) 



(13) 



- (14) 
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Eq. (15) cannot be solved explicitly, but an implicit solution can be derived instead. Note that Eq. (15) 
can be expressed as 



Integration of Eq. (16) gives 



Q'iO - Ql 

This equation can be integrated as follows: 
"1 dX 



dC 7 + 1 rfC 
202 



= 



7 + 



-(c^+i(C) -C7+') = 



(16) 



(17) 



C 



c 



'g? + ^(X7+i_c7+i) 

7 + 1 



12^ 
7 + 1 



- Ql 



c 



Q?(7 + l) 



dX 



I 202 
7 + 1 



- Ql 



202 

1 C'^+i 



(18) 



2'7+l' 7 + l'^7+i (7 + l)Q! 

^1 202 



2^i( ^,^-t;i 



\2'7+l' 7+l'^7+i (7 + l)Q^ 
^1 ^0^ 

where 2-^1(0, ^; c; a;) is the hypergeometric function. An analytical continuation formula is required to 
compute the hypergeometric function [22, 23]. Given that is a real number, it follows that 



Clm 



/II 1 
2Fi{ -, r;l ■ 



27 + 1' i + h + i)Qi 



-Im 



/II 1 
2F1I -, -■! 



202 
1 



2'7 + l'^ ' 7 + l'^7+i (7 + l)Qj 
^1 202 



(19) 



This solution holds for 202c7+V(7 + l)-Ql>0. If 

202^7+7(7 + i)-Q? = o 

it is obtained [21] 



1 



(I-7) 202 



7 + 1 



This is a special case where an explicit relation between C and arises. For 7 = 1 wc have 

(7(C) - cxp(-0C) 

If 

202C7+V(7+1)-Q?<O 



(20) 



(21) 
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then we can integrate Eq. (17) to give 



2Fi\ 



1 1 



2'7+l' 7+1' il + l)Ql 



202 



1 



1 



2'7+l'^ ' 7 + 1' (7 + l)Qi r^+i 



7 + 1 



(22) 



The profile of methanol coneentration can be computed from the above equation by using an iterative 
procedure. In order to calculate the anodic polarization curve, it is necessary another relation between the 
local anodic current density and the local methanol concentration at the CL. From Eq. (7) 



By invoking the chain rule: 



we obtain 



— = -g(C)=iJ(C)-Qo 

dC 



dJ _ dJ dC 



dJ 



.C7 



dC I 

noting that J(0) = when C = 1, an integration of the above equation gives: 

202 



I (7 + 1)1 
solving this second degree polynomial we obtain 



(1 -C+i) = 



■HO 



I- ^ l- 



202 



(1 - c-'+HO) 



(23) 

(24) 
(25) 

(26) 
(27) 



This equation relates the local anodic current density with the local methanol concentration at the CL. 
Also, Eq. (27) allows to express local anodic current density as a function of the anodic overpotential, 
kinetic and mass transfer parameters. This equation is valid if the fast proton transport assumption is 
fulfilled. Substitution of this equation into Eq. (23) gives 



Q(C) _ 

Oo 

At the CL/DL interface Eq. (17) becomes 



202 



(7 + 1)0^ 



(1- (77+1(0) 



7+1 



7 + 

It can be noted that Eq. (7) can be written as 

I = Qq — Qi 

Substitution of this equation into Eq. (29) and rearranging we obtain 



(7 + l)0g 



2(1 - C7+') 



1-1 



I \2 



(28) 



(29) 



(30) 



(31) 
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Substitution of this equation into Eq. (12) gives the relation (using the original variables): 



7 + 1 



7 + 1 



(32) 



2.3 Limiting cases 



For high cell current density regime where Cy^ 
expression arises 

lim (j)^ 

Ci->-0 

now it can be seen that cjp ^ oo as C^^^ 



■c\/n\ 



0, methanol flux crossover tends to zero and the following 



7 + 1 



RT 
aj- 



In 



7+1 



0, and the anodic overpotential 

/ 



nyiJ-kf-iC-y/ 



(33) 



(34) 



tends to infinity. This approximation avoids the calculation of the hypergeometric function and the iterative 
calculations for high fuel cell current densities. If Ci — > 0, Eq. (17) can be integrated to give Eq. (20). 
Therefore, Eq. (20) applies for high current densities and for the special case where 2(/>2C'7+V(7 + 1) - 0? = 
0. 



2.4 A simple model for methanol crossover 

In order to simulate dimensionless quantities, the methanol crossover by electro-osmotic drag is neglected. 
Thus, 

Qi - j^Ci = rCi (35) 

feci 



here 



-r)m 

k - -M. 

Om 



(36) 



where is the methanol diffusion coefficient in the membrane and 5^ is the diffusion layer thickness. This 
expression is used only for computational purpose; the ratio r was taken as 0.08 [12]. Electro-osmosis and 
pressure gradients must be included for realistic simulations. 



3 Discussion 

In order to illustrate the analytic method presented in this work, curves of the Thiele modulus ((/)^) for 
7 = 0.5 and 7 = 1 are shown in Fig. 2. As expected, the Thiele modulus tends to increase quadratically in 
the mass transfer region (i.e., (/) ^ 4) where the methanol crossover flux tends to zero. In this region, the ap- 
proximation (Eq. (34)) shown in Sec. 2.3 is applicable in this model, and easily allows to estimate the anodic 
polarization curve for high fuel cell current densities. For low and moderate cell currents (i.e., < </> < 2), 
Eq. (32) must be used. This Eq. gives the anodic overpotential as a function of methanol concentration 
at the CL/membrane interface and cell current density. Due to the fact that the hypergeometric function 
in Eqs. (19) and (22) cannot be inverted (except for some special cases), it is not possible to express the 
methanol concentration at the CL in terms of position. This implies that the methanol concentration at the 
CL/membrane interface cannot be expressed as an explicit function of the fuel cell current. Now, given that 
the anodic overpotential in Eq. (32) depends on the methanol concentration at the CL/membrane interface, 
it cannot be expressed as an exact explicit function of the fuel cell current uniquely. An alternative is to 
neglect the methanol concentration at the CL/membrane, suggesting that the Eq. (34) could be used to fit 
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experimental data. However, because Eqs. (34) and (32) do not correspond one to another, neglecting the 
methanol concentration at the CL / membrane interface could limit the applicability of Eq. (34) for fitting 
data. 

Fig. 3 shows plots of dimensionless methanol concentration at the CL/membrane interface vs. Thiele mod- 
ulus (0). It can be noticed that the decrease of dimensionless methanol concentration at the CL/membrane 
interface is more pronounced as the reaction order is reduced, giving a reduction of methanol crossover. Fig. 
3 can be interpreted as the diffusive part of the dimensionless methanol crossover flux. For a more detailed 
analysis of methanol crossover flux, electro-osmosis and pressure gradients must be added to the diffusive 
contribution of the methanol crossover flux [11, 13]. 

The dimensionless profiles of local anodic current density for various Thiele modulus are shown in Fig. 4, for 
moderate current densities, the profile shows a linear behavior where the methanol concentration in the CL 
is appreciable and nearly constant across the CL. This suggests that the proton generation is homogeneous 
in the CL. The local anodic current density becomes non-linear in the concentration overpotential region. 
In this case, the methanol concentration is negligible in the vicinity of the membrane, and the proton flux 
is generated mainly in the vicinity of the DL, and it becomes constant near the CL/membrane interface. 
Wang et al. [9] reported similar results for the dimensionless profiles of current density for the case of oxygen 
transport in a cathodic catalyst layer by assuming oxygen transport through spherical agglomerates and fast 
proton transport. 

The profiles of local methanol concentration at the CL are shown in Fig. 5, in the activation region there is 
a little concentration difference between C = and C = 1- In accordance with the results shown in Fig. 4, 
Fig. 5 show that in the mass transport limitation region, the methanol concentration at the CL/membrane 
interface and methanol crossover (Fig. 3) becomes negligible [11]. In this case Eq. (20) is suitable to calcu- 
late local methanol concentration at the CL profile. 

Finally, the Fig. 6 shows the profiles of dimensionless electrochemical reaction rate in the CL, these profiles 
can be related to the CL activity, for low current densities the activity is uniform across the CL, for high 
current densities, most of the activity of the CL is located near the DL/CL interface, and it controls the fuel 
cell performance, this confirms the characteristics of the Figs. 4 and 5. This analysis of the electrochemical 
reaction rate profiles is applicable for the reduction of oxygen at the cathode of a PEMFC [14] . 
Eq. (15) relates the methanol reaction and diffusion in the CL of the DMFC. The analytical solutions of this 
equation (Eqs. (19), (20) and (22)) are in general implicit functions (except Eq. (20)) of the local methanol 
concentration, the analysis of these equations is not easy, and it requires an iterative procedure to calculate 
the methanol profile in the CL. Once calculated the dimensionless methanol concentration profile, the local 
anodic current density can be computed by using Eq. (27). This Eq. gives a simple relationship between 
current density and the methanol concentration. It can be noted that Eq. (27) apparently shows that the 
profile of local anodic current density is non-linear. However, a Taylor expansion of Eq. (27) can show that 
the local anodic current density is a linear function of C for low current densities. Therefore, in the special 
case of low cell current densities, the profiles of local methanol concentration at the CL and local anodic 
current density can be expressed as linear functions of C. 

Fast proton transport {drj/dz = 0) was the key assumption to derive Eq. (19), if the electrical resistance of 
the CL is negligible, then the local anodic overpotential is nearly constant; this occurs when the CL is very 
thin and its pores are well covered by the SPE. Consequently, the analytical solutions presented in this work 
cannot be applied if the fast proton transport assumption is not fulfilled. For the case of oxygen reduction, 
the simulation presented by Bernardi and Vebrugge [14] shows a nearly constant overpotential profile in the 
cathodic catalyst layer, and Gurau et al. [15] take advantage of this fact to derive analytical solutions in 
order to simulate a porous electrode with 7=1. 

The analytic model for the anodic catalyst layer presented here is subjected to the limitations of other one- 
dimensional models (see [10, 16, 17]) that do not consider two-phase transport. For example, Jeng and Chen 
[10] modelled the anodic catalyst layer with Tafel Kinetics and two-phase transport was neglected. In their 
model, like ours, one-dimensional transport of methanol in the anodic catalyst layer was considered. They 
showed that their model accurately fits the anodic polarization curve for 2 M of methanol feed concentration. 
Nevertheless, further validation of their model was not reported. 
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Another model that considers spatial variations of methanol concentration in the anodic catalyst layer of 
a DMFC was reported by Garcia et. al [16]. But, in comparison with the model of Jeng and Chen [10], 
Garcia et. al [16] modelled the anodic catalyst layer of a DMFC using the kinetic expression of Meyers and 
Newmann [20], and it was assumed that the local anodic current density and the local anodic overpotential 
were independent of position. The simulation of the polarization curve reported by Garcia et. al [16] shows 
good agreement with the experimental data only for 0.05, 0.1, 0.2 and 0.5 M of methanol feed. They found 
that for methanol feed concentrations greater than 0.5 M it was necessary another set of fitted kinetic pa- 
rameters to model the polarization curve, and they suggest that this was caused by neglecting two-phase 
transport. In fact, it has been shown that models without gas transport do not agree well with experimental 
data for high fuel cell current densities [11. 13]. Therefore, in order to perform more accurate simulations, 
important phenomena like two-phase transport [11, 13], transport of species through the electrolyte [6, 18], 
energy transport [18], mass transport along the flow channels [19], and non-Tafcl kinetics must be included 
in a rigorous DMFC model. Although these phenomena were not considered in this work, the equations 
derived here can be used for rapid estimations, approximate quantitative analysis of the performance of the 
DMFC, and derivation of semi-empirical equations. 

4 Conclusions 

In this work, analytical solutions for the anode side of a DMFC performance simulation were presented. 
These solutions are applicable to the case of negligible CL proton resistance and one-dimensional transport. 
However, the equations obtained for the methanol concentration and local anodic current density profiles 
arc not explicit functions of the position in the CL. This implies that an iterative procedure must be used 
to obtain these profiles. 

A set of analytical equations that relate the local methanol concentration at the CL with dimensionless 
position was found. In general, these equations do not allow to find an expression of the local methanol con- 
centration at the CL as an explicit function of dimensionless position. Because of this, analytical expressions 
for local anodic current density and methanol flux at the CL as function of local methanol concentration at 
the CL were derived. Therefore, if the profile of local methanol concentration at the CL is known, then the 
local anodic current density and methanol flux at the CL profiles can be easily calculated by means of Eqs. 
(27) and (28) respectively. 

Thiele modulus was found to be a function only of mass transport quantities such as the methanol flux 
in the DL and the methanol concentration in the CL/membrane interface. This dependency is caused by 
the assumption of fast proton transport, if this assumption is not fulfilled, the expressions derived for the 
Thiele modulus are no longer applicable, and there is no analytical method available to simulate the cell 
performance; hence, a procedure based on numerical methods would be required instead. 
The expression derived in this work for the Thiele modulus allows the calculation of the anodic overpotential, 
and the profiles of methanol concentration and current density in the CL. Also, by means of the general 
expression for the Thiele modulus derived (Eq. (31)), the case of high fuel cell current densities was studied. 
A simple expression that relates the anodic overpotential with cell current density was obtained (34). This 
expression could be used to fit experimental polarization curves, but it might not work for low cell current 
densities. 
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A Notation 



C dimensionless methanol concentration as function of dimensionless position in CL 

Co dimensionless methanol concentration at the CL/DL interface 

Ci dimensionless methanol concentration at the CL/membrane interface 

methanol concentration in CL as fmiction of position (mol/cm"^) 

C^'''^' methanol concentration at the CL/DL interface (mol/cm"^) 

C^^™ methanol concentration at the CL/membrane interface (mol/cm'^) 

CJ^* refer ence methanol concentration (mol/cm'^) 

Dm diffusion coefficient of methanol in water (cm^/s) 

-^M*^ effective diffusion coefhcient of methanol in CL (cm^/s) 

T>y^ diffusion coefficient of methanol in the membrane (cm^/s) 

Faraday constant (C/mol) 

/ cell current (A/cm^) 

I dimensionless cell current 

j local anodic current density (A/cm^) 

reference exchange current density (A/cm^) 

J dimensionless local current density 

fcci mass transfer coefficient of methanol in the CL (cm/s) 

fcm mass transfer coefficient of methanol in the membrane (cm/s) 

riM number of electrons transferred per methanol molecule 

Q dimensionless methanol flux as function of dimensionless position 

Qo dimensionless methanol concentration at the CL/DL interface 

Q\ dimensionless methanol concentration at the CL/membrane interface 

A^;^ local methanol flux in CL (mol/cm^ s) 

Ny^ methanol flux in DL (mol/cm^ s) 

methanol flux in M (mol/cm^ s) 

R ideal gas constant (J/molK) 

S CL specific area (cm^/cm^) 

T cell temperature (K) 

V convective velocity 

z position in CL 

a charge transfer coefficient 

7 reaction order of methanol electro-oxidation 

6ci thickness of CL (cm) 

5m thickness of m (cm) 

Cci CL void fraction 

f] anodic overpotential (V) 

0^ Thiele modulus 

^ dimensionless position in CL. 
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Figures: 

Figure 1. Schematic representation of the anode side of a DMFC. 

Figure 2. Plots of Thiele modulus vs. dimensionless current for 7 = 0.5 and 7 = 1. 

Figure 3. Plots of dimensionless methanol concentration at the CL/membrane vs Thiele modulus for 7 — 0.5 
and 7 = 1. 

Figure 4. Dimensionless methanol concentration profiles for various Thiele modulus; 7 ~ 0.5. 
Figure 5. Dimensionless local anodic current density profiles for various Thiele modulus; 7 = 0.5. 
Figure 6. Electrochemical reaction rate profiles for various Thiele modulus; 7 = 0.5. 
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